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Abstract 



2 OUR STRATEGY 



A practical strategy is presented and successfully imple- 
mented to determine form factors in the time-like but un- 
physical (below threshold) region using dispersion rela- 
tions, in a model independent way without any bias towards 
expected resonances. Space and time-like data have been 
employed along with a regularization scheme to unfold and 
solve the integral equations. Remarkably, resonance struc- 
tures with peaks for the p{770), p'(1600) and a structure 
near the NN threshold are automatically generated. The $ 
peak is invisible thus refuting suggestions about any size- 
able ss content in the nucleon. 

1 THE "IN PRINCIPLE" METHOD 

Consider a dispersion relation (DR) for a (generic) normal- 
ized nucleon form factor (FF) G{s) with a subtraction at 
s = and t <0 
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For s > So (the lowest mesonic threshold, 4to^ for even 
and 9m^ for odd G parity NN channels), the phase 6{s) 
of the FF defined as G{s) = \G{s)\e'^^'''> is given by the 
principal value integral 



Sis) 



s^/s'- 



■Pr 



ln\G{s')\ds' 
s'{s' - s)y/P^ 



(2) 



So in principle the method is very simple: (i) use space-like 
data on the left side of Eq.(l) and solve the integral equa- 
tion to find In |G(s)| for s > Sq. (ii) Having determined the 
modulus hi |G(s)|, use Eq.(2) to compute the phase S{s). 

The above procedure has been unsuccessful in the past as 
it is an "ill-posed" mathematical problem[|l], |[]. The result 
depends on the input data in an unstable way and an im- 
possible accuracy is needed before one arrives to a stable 
unique solution. 



A successful method||3| has been developed by splitting the 
time-like region into 2 parts: 

(i) Region I : is the unknown unphysical region [sg, 4m^] 
for which the FF is to be determined. 

(ii) Region II : consists of the physically accessible time- 
like region s > Amj^, for which data exist and quite accu- 
rate asymptotic estimates are available. 

With the above breakup, the unknown part of the inte- 
gral equation is reduced to the (small) region I, which is 
amenable to a finite matrix analysis with some technical 
refinements described below in brief (the details of the de- 
veloped procedure may be found in references and [^). 

An integral equation of the first kind, linear in the un- 
known In |G|, can be derived 

InG(t) - m ^ f^-\^n\Gis)\d^^ 



where 
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is a "known" quantity since it can be calculated directly 
from experimental data in the time-like region with some 
recipe to extrapolate them to very high t values. 

To avoid instabilities usually met in solving first kind 
integrals such as Eq.(3), we impose a regularization scheme 
with smoothness: 

• lit) is calculated using time-like data through a ratio- 
nal, smooth function with the expected asymptotic behav- 
ior. The subtraction at s = helps in diminishing the im- 
pact which the asymptotic behavior has on the results. 

• There is a steep spike near the NN threshold. To avoid 
any ensuing instabilities, the upper limit of the unphys- 
ical region has been raised to S2 = 4m^ + A where 
A « 0.5GeV^ and continuity is imposed _there. A new 
DR is constructed for the region (4m^, S2) 



Figure 1 ; Pion space-like FF computed via DR with- 
out (grid band) and with (light band) subtraction, 
compared with pion FF space-like data and square 
route of proton space-like FF data fit (dark band). 
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Figure 3: Phases of proton (light band) and pion (dark 
band) FF computed via DR. 



Figure 2: Pion FF computed via DR from time-Uke 
data. 
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Figure 4: Phase of pion FF (dark band) compared 
with the half of nucleon iso vector FF (light band). 



• Our regularization consists in requiring the local cur- 
vature of the FF in the unphysical region, i?2 = 
Is^ (^^-^5^)^rfs to be limited. Instead of the second 



derivative of In jCj, as is standard|[lj we employ the sec- 
ond derivative of \G\ for this purpose. The reason is that 
fluctuations in \G\ are important only when \G\ is large, 
while In |G| fluctuations would be large also when |G| is 
small. 

• Eq.(3) is then linearized by transforming the integrals into 
sums over M — 50 suitable sub intervals in s, with their 
widths increasing with s. This introduces further smooth- 
ness, by effectively integrating over any structure with a 
narrower half width. 



• The minimize the integral 
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where — ln|G[(sj+i +Sj)/2]| is calculated in the mid- 
dle of the 2 th sub interval, ti, with j = 1, • • • , L, corre- 
spond to the experimental points available in he space-like 
region. 

• Finally, we have 



T^total = Ro + T^i?2 + C* 
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Figure 5: Proton magnetic FF computed via DR. 
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Figure 7: Imaginary part of the nucleon magnetic 
isovector FF computed via DR compared with expec- 
tation from unitarity relation. 
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Figure 6: Neutron magnetic FF computed via DR. 
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Figure 8: Imaginary part of the nucleon isoscalar FF 
computed via DR. 



The "dumping parameter" t has to be chosen by trial and 
error: if it is set too large, it will not respond to sharp struc- 
tures, while unstable solutions will result if it is set too low. 
• The uncertainties in the solution of Eqs.(2) and (3), due 
to experimental errors, were estimated by simulating new 
space- and time-like data according to the quoted errors and 
then solving the DR for each simulated set. 

3 TEST OF THE REGULARIZATION 
METHOD 

To test the entire procedure and also to get a suitable range 
for the parameter t, we computed the space-like pion FF 
using time-like data. In the time-like region, this FF is 
known up to the J/^P mass and at higher Q^, it was extrap- 
olated using first order QCD[^. In Fig.(l), comparison is 
made with the measured (low Q^) space-like data. (Higher 



space like data points are through extrapolations from 
pion electroproduction data and thus may have systematic 
errors). We have also made other tests[3,4] obtaining good 
agreement with the p peak, the p width and also a dip at 
1.6 GeV^ for r w m.^ Fig. (2). The phase of the pion FF 
approaches just above 2 GeV to its expected asymptotic 
value of 180 degrees Fig. (3,4). 

4 RESULTS FOR THE NUCLEON FF 

Below we summarize some salient features of our findings: 

• For the first time, resonant structures have been generated 
from "smooth" inputs Fig. (5, 6,7, 8). The method is stable 
and reliable. 

• The combined {p + lo) peaks and the p'(1600) are gen- 
erated at the right mass. However, the p peak is much 
broader Earlier analyses!^ had also found a similar dis- 



crepancy. 

• No $ signal is visible thus signaling a very small ss con- 
tent in the nucleon. 

• Phases for the nucleon are consistent with expectations: 
Sn 360° within the error bands Fig. (3). 

• There is an interference pattern near threshold (M « 
1.88 GeV) which may be related to baryonium Fig. (5, 6). 

• ImG^j^' changes sign once and ImG^j' appears to change 
sign twice and various superconvergence sum rules are all 
obeyed by our FF's Fig. (7,8). By way of comparison, our 
analysis strongly indicates that Imi^^r does not change sign. 
Neglecting logaritmic factors (which we can not resolve), 
this would suggest (for power law behavior) that 

|F^(s)| |s|-i+' as \s\ oo. (e > 0). 

5 CONCLUSIONS 

Nucleon time-like magnetic FF have been obtained in an al- 
most model independent way by means of DR for In G{q^), 
using a regularization method in conjunction with space 
and time-like data. Resonances have been found consis- 
tent with the p(770) and p'(1600) masses. However, a 
very large p width is obtained - as in previous DR anal- 
yses. Further work is in progress to understand the sources 
of discrepancies as well as the relationship of our results 
with other DR analyses |^]. 

Other applications of this strategy have been discussed at 
this conference by Y. Srivastava (see contribution T19). 
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